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EXPLICIT  SOLUTION  OF  DIFFERENCE  EQUATION  FOR  THE 
WAVENUMBER  RESPONSE  OF  FLUID-LOADED  STIFFENED  PLATE 


INTRODUCTION* 


Acoustic  radiation  from  the  upper  surface  of  a  fluid-loaded 
stiffened  plate  has  been  examined  by  Cray  [1].  The  formulation 
developed  there  produces  an  implicit  algebraic  equation  for  the 
Fourier  transformed  plate  wavenumber  response.  This  response 
must  be  manipulated  into  an  explicit  form  in  order  to  obtain  a 
solution  for  the  plate's  structural  response  and  to  obtain  the 
near-  and  far-field  generated  acoustic  response.  The  following 
is  a  brief  synopsis  of  the  formulation  for  the  acoustic  radiation 
from  a  stiffened  plate. 

The  infinite  plate  investigated  in  [1]  was  configured  to  have 
two  infinite  sets  of  attached  rib  stiffeners.  The  stiffeners 
composing  a  given  set  were  identical  and  were  spaced  periodically 
with  distance  L.  One  set  of  stiffeners,  however,  was  shifted  by 
an  amount  A  from  the  other  set.  In  this  manner,  portions  of  the 
plate  were  configured  with  repeating  sections  having  nonperiodic 
rib  spacing. 

The  governing  equation  of  motion  for  the  surface  displacement 
of  the  fluid-loaded  isotropic  plate  for  an  applied  line-force, 
with  Q  sets  of  attached  rib  stiffeners,  is  given  by 


—4--^  -  m  w(x)  «  P^(x)  -  P,(x,0)  -  ^  P„(x 


(i: 


dx" 


q-1 


*The  introduction  was  contributed  by  Dr.  Benjamin  A.  Cray. 
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where  D  is  the  rigidity  of  the  plate^  w(x)  is  the  transverse 
plate  displacement/  m  is  the  mass  per  unit  area  of  the  plate, 
fa)  is  the  applied  excitation  frequency,  Pg(x)  is  the  external 
pressure  due  to  the  applied  line-force,  P,(x,0)  is  the  acoustic 

a 

pressure  acting  on  the  upper  surface  of  the  plate,  and  {Pq(x)l 
are  the  total  pressures  exerted  by  each  set  of  attached  rib 
stiffeners . 

Rotary  inertia  and  shear  deformation  effects  within  the  plate 
and  rib  stiffeners  were  neglected.  The  stiffeners  exerted 
reactive  forces  upon  the  plate  but  applied  no  angular  moments. 

Equation  (1)  was  transformed,  term  by  teirm,  using  exponential 
Fourier  transforms.  The  spatial  transform  variable,  k,  has  the 
physical  significance  of  wavenumber.  It  was  assumed  that  the 
transforms  are  well  defined  and  exist  over  the  entire  domain  of 
integration. 

Upon  transforming  (1)  into  the  wavenumber  domain,  the 
implicit  form  of  the  wavenumber  response,  for  the  fluid-loaded 
stiffened  plate,  can  be  written  as 


00 

W(k)  =  F(k)  -  Y(k)  YIZ  W(k  +  nkj^)  x 

ns-oo 

X  (bj  +  B2  explink^hj  +  .  .  •  +  b^  exp (inkLhQ_  J  )  , 


where  W(k)  is  the  (transformed)  wavenumber  response.  The 
quantities  on  the  right-hand  side  of  (2)  are  defined  as 

F  exp(-ikx  )  j 

-  sW) — 


(2) 
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where 


S(k) 


-  ''b) 

b(b^  -  bj) 


for  |k|  <  k 

for  |k|  >  k 
'  '  o 


(4) 


2  k 

where  kj^  =  (m  u  /D)  is  the  free  in-vacuum  plate  wavenumber, 
k^  =  2n/L  is  the  wavenumber  defined  by  the  periodic  spacing,  and 
Bq  =  m^  (ii  /L  is  the  dynamic  structural  mass  of  the  q-th  rib 
stiffener  set.  Also,  is  the  q-th  offset  distance. 

Notice  in  (2)  that  the  wavenumber  or  spectral  response  W(k) 
appears  implicitly.  The  difficulty  now  lies  in  determining  an 
explicit  expression  for  the  plate's  wavenumber  response  W(k). 

That  is,  it  is  necessary  to  manipulate  (2)  such  that  the 
summation  which  contains  the  shifted  wavenumber  responses 
lW(k  +  nkj^)  )  may  be  rewritten  in  terms  of  known  quantities. 

In  [1],  Cray  obtained  an  explicit  expression  for  W(k)  for  the 
case  of  a  single  offset,  that  is,  Q  =  2  in  (2).  This  corresponds 
to  two  sets  of  rib  stiffeners,  where  one  set  is  offset  from  the 
other.  The  following  development  is  a  significant  extension  of 
the  configuration  of  two  sets  of  rib  stiffeners  to  the  case  of  an 
arbitrary  number  of  rib  stiffener  sets,  each  with  different 
offsets  {Aq},  and  to  arbitrary  excitation  functions. 
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DESCRIPTION  AND  SOLUTION  OF  DIFFERENCE  EQUATION 

In  the  following  mathematical  development,  k  and  t  are 
arbitrary  real  variables,  while  m,n,p,q  are  integers.  Also, 
sximmations  without  limits  are  from  -®  to  +®.  We  are  interested 
in  finding  the  solution  W(k)  to  the  difference  equation 

W(k)  =  F(k)  -  Y(k)  Y.  \  W(k  +  cn)  ,  (5) 

n 

where  F(k)  and  Y(k)  are  given  known  functions  of  k,  and  sequence 
{A^l  is  also  known.  For  e  =  0,  the  solution  to  (5)  is  immediate; 
hence,  e  is  nonzero  and  fixed  in  the  following. 

In  fact,  we  are  interested  in  the  particular  case  where  A^^  is 
a  finite  sum  of  Q  complex  exponentials  in  n  (compare  with  the 
second  line  of  (2)): 

A„  =  >  '  exp(ina  )  for  all  n  ,  (6) 

where  {otg),  1  1  q  1  Q,  are  arbitrary  (complex)  distinct 
constants,  and  {B  }  are  arbitrary  complex  constants.  The  (B  ) 
and  (0^1  are  given.  Furthermore,  since 

exp(inaq)  =  exp(in[aq  -  p2n])  for  all  p  ,  (7) 

there  is  no  loss  of  generality  in  assuming  that  |Re(ag)|  S  n  for 
all  1  ^  q  1  Q.  If  Im(aq)  is  nonzero  for  any  q,  then  [a^^I  cannot 
remain  bounded  at  both  limits  of  n  «  ±®.  The  earlier  case  solved 
in  [1]  corresponded  to  Q  »  2  and  the  special  value  of  ®  0. 
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DEFINITION  AND  PROPERTIES  OF  AUXILIARY  FUNCTIONS 
Define  functions 

W  (k)  =  X]  W(k  +  en)  exp(ina  )  for  1  i  q  £  Q  .  (8) 

^  n  ^ 

(The  additional  dependence  of  Wq(k)  on  e  is  suppressed 
notationally • )  It  then  follows  that 

W  (k  +  cm)  =  5IW(k  +  Ein+  En)  exp(ina  )  = 

^  n  ^ 


=  W(k  +  Ep)  exp(i{p-m)a  )  =  exp(-ima  )  W^Ck)  ,  (9) 

P  q  q  -q 

where  we  let  p  =  m  +  n  and  used  (8).  This  is  the  key  relation 
regarding  the  functions  Wq(k)  defined  in  (8),  namely 

(k  +  Em)  =  exp(-ima  )  W  (k)  for  1  £  q  £  Q  .  (10) 

with  the  aid  of  (6)  and  (8),  the  summation  on  n  in  (5)  can 
now  be  manipulated  into  the  form 

W(k  +  En)  =  X]  +  En}  B  exp(ina  )  =  X^  B  W^(k)  . 
n  "  n  q=l  ^  q=l  ^ 

(11) 

This  enables  us  to  express  (5)  in  the  alternative  form 
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SOLUTION  OF  EQUATION  (12) 

At  this  point,  in  analogy  to  (8),  it  is  convenient  to  define 
two  additional  functions,  namely, 

F  (k)  =  Y.  exp(ina  )  for  1  i  q  i  Q  ,  (13) 

^  n  ^ 

Xpq(*^)  =  E  Y(k  +  en)  exp[in(ap  -  Og)]  for  1  i  p,q  i  Q  .  (14) 

Now,  replace  k  by  k  +  cm  in  (12)  and  use  (10),  obtaining 
W(k  +  cm)  =  F(k  +  cm)  -  y(k  +  cm)  >  ;  h  exp(-ima  )  W  (k)  .  (15) 

q  r-i  q/  _qv 

When  this  equation  is  multiplied  by  exp(imap)  and  summed  over  all 
m,  there  follows,  by  use  of  (8),  (13),  and  (14), 

Wp{k)  =  Fp(k)  -  Ypg(k)  Bg  Wg(k)  for  1  £  p  i  Q  .  (16) 

This  relation  constitutes  Q  linear  equations  (at  each  k)  in  the  Q 
unknown  functions  Wq(k),  1  <  q  £  Q.  The  remaining  quantities  in 
(16)  can  be  obtained  from  (6),  (13),  and  (14). 

When  (16)  is  solved  for  all  the  {Wq(k)l,  then  (12)  directly 
yields  the  original  quantity  of  interest,  namely  W(k)  in  (5),  for 
that  particular  value  of  k.  Although  this  procedure  appears  to 
require  the  solution  of  a  new  set  of  Q  simultaneous  linear 
equations  for  each  k  of  interest,  there  is  a  shortcut  that  might 
be  useful  in  some  cases.  Namely,  reference  to  (15)  reveals  that 
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W(k  +  em)  can  now  be  calculated  easily  for  nonzero  m,  for  that  k, 
once  F(k  +  cm)  and  Y(k  +  em)  have  been  computed,  without  the  need 
for  another  solution  set.  The  same  solutions  {Wq(k)}  are  used  on 
the  right-hand  side  of  (15),  regardless  of  the  value  of  m  under 
consideration;  only  the  finite  summation  in  (15)  need  be  redone. 
Thus,  the  solution  W  to  (5)  at  arguments  k,  k  ±  e,  k  ±  2e,... 
can  be  found  from  the  solution  of  one  set  of  Q  simultaneous 
equations,  (16). 

For  the  special  case  of  Q  =  1,  (16)  can  be  immediately  solved 
to  yield  [ 1 ] 


P  H  F(k  +  En)  exp(inaj) 

=  1  Y  -  • 

1  ^  ®1  -11^*^^  1  +  Bj  Z  Y(k  +  En) 

n 


Then  (12)  yields  the  explicit  result  (for  Q  =  1) 


F(k  +  En)  exp(inaj) 

W(k)  =  F(k)  -  B,  Y(k)  -  .  (18) 

1  +  XI  +  En) 

n 

In  this  form,  it  is  possible  to  directly  evaluate  W(k)  at  any  k 
of  interest  merely  by  computing  the  terms  encountered  on  the 
right-hand  side  of  (18).  However,  shortcut  (15)  is  still  better 
for  the  particular  arguments  (k  +  Em). 
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For  Q  =  2,  when  the  solutions  of  (16)  are  substituted  into 
(12),  the  function  W(k)  is  given  by  the  explicit  expression 


W(k)  =  F(k)  -  Y(k)  , 


(19) 


where 


N(k)  =  Bj  Fj(k)  +  B2  F2(k)  +  Bj  B2  X 

^  [ll(k)  X22<^)  X21<^)  "  ^^2^’^)  ^12^’^^] 


(20) 


and 


D{k)  =  1  +  Bj  Y^j(k)  +  B2  Y22(k)  + 

+  Bj  B2  [Yjj(k)  Y22(k)  -  Yj2(k)  Y2i(k)] 


(21) 


A  program  for  the  evaluation  of  (19)  -  (21),  for  the  case  where 
Y(k)  is  real,  is  presented  in  the  appendix;  also,  B^,  B2/  and 
02  are  real.  Then  D(k)  is  real,  but  N(k),  F(k),  and  W(k)  are 
complex.  Also,  y22(^)  ~  — 11^^^  real,  while  Fj(k),  F2(k),  and 
Y2i(k)*  =  Yj2(k)  are  complex.  The  program  uses  these  properties. 

In  general,  for  any  Q,  if  (16)  is  solved  analytically  for  the 
set  (Wq(k)}  in  terms  of  {Bg},  (Fq(k)),  and  iYpq(k)l,  these 
results  can  be  substituted  or  utilized  in  (12)  to  get  an  explicit 
expression  for  W(k)  that  is  valid  for  all  k.  For  large  Q,  this 
will  be  impractical;  instead,  numerical  solution  of  the  Q  simul¬ 
taneous  equations,  (16),  will  be  required  at  each  k  of  interest. 
The  only  exception  is  to  use  (15)  for  the  particular  arguments 
(k  +  cm}.  (All  results  above  actually  hold  for  k  and  c  complex.) 

9/10 
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SOLUTION  FOR  GENERAL  SEQUENCE  {A^j 

In  this  section,  we  no  longer  restrict  sequence  {A^j  to  have 
exponential  form  (6).  Rather,  for  bounded  given  sequence 
define  function 

B(u)  =  ^  exp(-inu)  for  |u|  <  n  ,  (22) 

where  u  is  real.  Then,  the  sequence  values  can  be  found  from 

n 

A^  *  J  du  B(u)  exp(inu)  for  all  n  .  (23) 

-n 

The  fastest  variation  with  n  that  (23)  allows  for  A^  is 

exp(±inn);  however,  this  is  no  loss  of  generality,  as  seen  by 

reference  to  (7).  The  case  considered  earlier  in  (6)  corresponds 

to  B(u)  being  composed  of  a  set  of  Q  impulses  of  area  B^  located 

at  u  =  a„,  with  |a„|  In,  when  a  is  real, 
q  '  Q  Q 

Now,  define  function 

W(u,k)  =  W(k  +  en)  exp(inu)  for  |u|  S  n  .  (24) 

n 

(Again,  the  notational  dependence  on  e  is  suppressed.)  The  key 
relation  that  W(u,k)  satisfies  is 

W(u,k  +  em)  *  W(k  +  cm  +  cn)  exp(inu)  = 
n 

“  51  W(k  +  ep)  exp(i(p-m)u)  «  exp(-imu)  W(u,k)  ,  (25) 

P 

where  we  let  p  *  m  +  n  and  used  (24).  This  relation. 
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W(u,k  +  era)  =  exp(-iinu)  W(u,k)  ,  (26) 

is  the  analogue  of  (10)  earlier. 

The  summation  in  (5)  can  now  be  expressed  as 

n 

^AW(k+  en)  =  5Zw(k+  cn)  [  duB(u)  exp(inu)  = 

n 

=  J  du  B(u)  W(u,k)  ,  (27) 

-n 

where  we  used  (23)  and  (24).  Therefore,  the  equation  of 
interest,  (5),  now  takes  the  form 

n 

W(k)  =  F(k)  -  Y(k)  J  du  B(u)  W(u,k)  .  (28) 

-n 

We  now  define,  in  analogy  to  (13)  and  (14),  the  two  functions 

F(u,k)  =  Y1  exp(inu)  for  |u|  i  n  ,  (29) 

n 

Y(v-u,k)  =  exp(in(v-u))  for  |u|,|v|  £  n  .  (30) 

n 


When  we  replace  k  by  k  +  cm  in  (28) ,  and  use  (26),  there 
follows 


n 

W(k  +  cm)  -  F(k  +  cm)  -  Y(k  +  cm)  J  du  B(u)  W(u,k)  exp(-imu)  . 

■"  (31) 

We  now  multiply  (31)  by  exp(imv),  v  real,  and  sum  over  all  m,  to 
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obtain 

n 

W(v,k)  «  F(v,k)  -  J  du  Y(v-u,k)  B(u)  W(u,k)  for  |v|  n  ,  (32) 

-n 

where  we  used  (24),  (29),  and  (30).  This  is  a  linear  integral 
equation,  with  kernel  y(v-u,k)  B(u),  for  unknown  W(u,k);  that  is, 
for  each  k  of  interest,  a  new  linear  integral  equation  must  be 
solved.  Then,  W(u,k)  will  be  known  for  |u|  i  n  at  that 
particular  k. 

Once  W(u,k)  is  known,  (28)  gives  desired  solution  W(k)  at 
that  particular  k  value,  since  F(k),  y(k),  and  B(u)  are  known 
functions.  If  a  large  number  of  k  values  are  of  interest,  (32) 
can  involve  a  great  deal  of  computational  effort;  however, 
arbitrary  bounded  sequence  {A^j  is  now  allowed  in  (5). 
Furthermore,  (31)  then  yields  the  solution  for  W  at  arguments 
k,  k  ±  E,  k  ±  2e,...  without  the  need  to  solve  another  integral 
equation. 

The  linear  integral  equation  in  (32)  could  be  solved  by 
recursion;  for  example,  the  p-th  order  approximation  is  given  by 


n 

Wp(v,k)  S  F(v,k)  -  J  du  Y(v-u,k)  B(u)  Wp_i(u,k)  for  |v|  £  n  . 

(33) 

A  possible  starting  value  for  this  recursion  is  WQ(u,k)  =  F(u,k). 


13/14 

Reverse  Blank 


TR  10015 


SUMMARY 

When  the  governing  differential  equation  of  motion  for  the 
surface  displacement  of  the  fluid-loaded  isotropic  plate  for  an 
applied  line-force,  with  Q  sets  of  attached  rib  stiffeners,  is 
transformed  into  the  wavenumber  domain,  a  difference  equation  is 
encountered.  Solution  of  this  latter  equation  for  exponential 
excitation,  for  a  particular  value  of  wavenumber  k,  is 
accomplished  through  the  definition  of  auxiliary  functions 
involving  sums  of  displaced  wavenumber  responses  with  exponential 
factors.  Finally,  solution  of  a  simultaneous  set  of  Q  linear 
equations  completes  the  required  calculations.  For  the  most 
general  excitation,  the  set  of  equations  becomes  infinite,  that 
is,  a  linear  integral  equation  must  be  solved. 
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APPENDIX.  PROGRAM  FOR  (19)  -  (21) 


In  the  following  program  for  solution  W(k)  given  by 
(19)  -  (21),  specified  function  Y(k)  is  presumed  real;  however, 
F(k)  is  allowed  to  be  complex.  The  following  example  corresponds 
to 

F(k)  *  exp (-.71  k^)  +  i  exp(-.91  k^  -  k)  , 
y(k)  =  exp(-.63  k^  +  k)  . 


le 

K*.7  !  k,  uavenumber 

20 

E= . 4  !  eps i 1  on 

30 

Ll=.l  1  alphal 

40 

L2^.2  1  a1pha2 

50 

Bl  =  l .  1 

60 

B2*l .  3 

70 

CALL  Fqk (k,E,Ll ,L2,Flr,Fl i  ,F2r,F2i  ) 

!  (13) 

80 

CALL  Ypqk (K,E,L1,L2, Yll,yi2r,yi2i) 

!  (14) 

90 

B«B1*B2 

100 

YaYl 1-Y12r 

110 

D*Y1 1*Y1 l-Y12r*Y12r-Y12 j*Y12i 

120 

D=1.+<B1+B2)*Y11+B«D 

!  (21 ) 

130 

Nr»<Flr-t-F2r  )*Y-<F1  i  -F2i  )*Y12i 

140 

Nr=Bl*Flr+E2*F2r+B*Nr 

!  (20) 

150 

Ni *<F1 i +F2i )*Y+<Flr-F2r)*Y12i 

160 

Ni=Bl*Fli+B2*F2i+B*Ni 

!  (20) 

170 

PRINT 

180 

PRINT  "DENOMINfiTOR  =  ”;0 

190 

Y»FNY<KJ/D 

200 

CALL  F<K,Fr,Fi ) 

210 

Wr=Fr-Y*Nr 

220 

W1 *Fi -Y»N1 

!  (19) 

230 

PRINT  Wr,Wi 

240 

END 

250 

1 
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260 

SUB  Fqk<K,E,Ll, 

L2,Flr,Fli ,F2r,F2i  ) 

270 

DOUBLE  N 

!  INTEGER 

280 

Flr=Fl i =F2r*F2i 

S0  . 

290 

FOR  N«-40  TO  40 

300 

CftLL  FCK+E*N,Fr 

,Fi> 

310 

IF  flBS<N)=40  THEN  PRINT  ftBSCFr )+ftBS<Fi > 

320 

C»C0S<N*L1 ) 

330 

S«SIN(N*L1 ) 

340 

Flr*Flr  +  Fr*C-Fi'»S 

350 

Fli^Fli+Fr^S+Fi+C 

360 

C*C0S<N*L2) 

370 

S»SIN(N*L2) 

380 

F2r=F2r+Fr*C-Fi *8 

390 

F2i =F2i +Fr*S>F 1 *C 

400 

NEXT  N 

410 

SUBEND 

420 

1 

430 

SUB  YpqkCK,E,Ll 

,L2,Yll,V12r,Y12i) 

440 

DOUBLE  N 

•  INTEGER 

450 

Y1 l=Y12r=Y12i =0 

• 

460 

FOR  N=-40  TO  40 

470 

T»FNY<K>E*N) 

480 

IF  flBSCN)=40  THEN  PRINT  T; 

490 

ft*N*<Ll-L2) 

500 

Y1 1=Y1  MT 

510 

Y12r=Y12r+T*C0SCfl) 

520 

Y12i»Yl2MT*SIN<fl) 

530 

NEXT  N 

540 

SUBEND 

550 

! 

560 

SUB  F<K,Fr,Fi ) 

570 

Fr*F i »0 . 

580 

fl».  7MK*K 

590 

IF  fl>100.  THEN 

610 

600 

Fr»EXP<-H) 

610 

fi*.91*K«KfK 

620 

IF  fl)100.  THEN 

640 

630 

Fi -EXPC-fi) 

640 

SUBEND 

650 

1 

660 

DEF  FNYCK) 

670 

fi«. 63*K*K-K 

680 

IF  fl>100.  THEN 

RETURN  0. 

690 

RETURN  EXPC-fl) 

700 

FNEND 

!  <  1  3  > 


'  CM) 


!  FCk) 


!  real  YCk) 
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